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Abstract In Chandrayaan I Mission and Mars Orbiter Mission, ISRO put the space- 
craft from Earth Parking Orbit to the trans-lunar/trans-planetary injection point 
through a series of orbit-raising manoeuvres around Earth. The start time and dura- 
tion of each burn, along with the size of the intermediate orbits and the burn sequence, 
must be designed, such that at the end of manoeuvres, the trans-lunar/trans-planetary 
injection epoch is achieved. This study aims to design the optimal multiple finite burn 
orbit-raising manoeuvres for trans-lunar and trans-planetary manoeuvres, using dif- 
ferential evolution as the optimization technique. The reduction in the computational 
time is achieved through multi-threading. 
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1 Introduction 


ISRO has conducted several space missions, notably, Chandrayaan I Mission and 
Mars Orbiter Mission. The spacecraft which carried the experimental payloads for 
the two missions had to orbit around the Moon/Mars in Lunar/Planet Parking Orbit. 
There are a series of manoeuvres before the spacecraft reaches that orbit. The space- 
craft is first injected by the launch vehicle into Earth Parking Orbit (EPO). The 
apogee/energy of the orbit is gradually increased until the spacecraft is injected 
into the trans-lunar/trans-mars trajectory. The spacecraft is then captured in an 
orbit around Moon/Mars through another sequence of manoeuvres. The orbit-raising 
manoeuvres are executed through finite burns for a specified duration around perigee. 
The burns are directed using velocity steering, in which the thrust direction is contin- 
uously aligned in the direction of the velocity of the spacecraft. It is well known that 
finite burn results in loss in the propellant mass. The finite burn loss can be reduced 
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by executing the burn when the spacecraft is close to the perigee/apogee of the orbit, 
which makes the gravitational force perpendicular to the thrust direction [1]. 


2 Problem Formulation 


The orbit-raising manoeuvres by ISRO are executed using a 440 N liquid thruster. 
The velocity impulse to be imparted is of the order of 1 km/s. This results in large 
burn duration and hence, large finite burn losses for a single burn. If the velocity 
impulse is imparted in multiple burns near the perigee, it will reduce the finite burn 
loss. For Chandrayaan I Mission, energy requirements of up to 160 m/s could be 
saved, which translates to around 40 kg of payload for the Chandrayaan I spacecraft. 
Also, these finite burn manoeuvres must be carried out such that there is no unplanned 
apsis line rotation. Further, there can be additional constraints, which are specific to 
different missions. For example, Chandrayaan I Mission had an additional constraint 
that the spacecraft should be visible from any one of the ISRO ground stations during 
each burn. Designing the multiple burns that save propellant mass and satisfy the 
constraints is an important task. The burn sequence (orbit numbers), intermediate 
orbit sizes, burn start time and the duration of burn must be chosen to satisfy the 
constraints and to make the spacecraft at perigee at the trans-lunar/planetary epoch. 


3 Solution Procedure 


This study devises an optimal strategy that consists of sequence of the orbit-raising 
manoeuvres which makes the spacecraft to be at the trans-lunar/trans-Mars injec- 
tion point at the specified injection Epoch. This study formulates the problem as 
an optimization problem with unknown design parameters as the burn sequence, 
intermediate orbit size, the start time and duration of each burn. Differential evolu- 
tion (DE), an optimization technique, is used to determine the unknowns because 
of its effectiveness to find a global optimal solution. Differential evolution does 
not require an initial guess but finds solution from a specified set of bounds on the 
unknown design variables [2]. Because the problem involves numerical integration 
of equations of motion, including perturbations, the process of finding the global 
optimal solution under finite burn requires enormous computational time [3, 4]. 

To reduce the computational time, a novel scheme has been devised. In this 
scheme, initial estimates of some of the unknowns are generated first by assum- 
ing impulsive burn in the optimization process, and then, these initial estimates are 
used in the optimization process with finite burn to generate the refined estimates. 
The optimization process assuming impulsive burn generates initial estimates for two 
unknowns (i) the sequence of burns and (ii) intermediate orbit sizes. The initial esti- 
mates of intermediate orbit sizes are used to generate initial estimates of the start time 
and duration of each burn assuming symmetric burn around perigee. The objective 
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function for the optimal impulsive burn strategy is set as minimizing the deviations 
in the trans-lunar injection (TLI) epoch and the final target orbit elements. These 
estimates are used in the optimization process where the orbit raising is modelled 
with finite burn. The intermediate orbit sizes obtained in this optimization process 
are used as target intermediate orbits for the finite burn optimization. It is to be noted 
that the intermediate orbit sizes may undergo changes in finite burn strategy, due to 
visibility constraint. 

The initial estimates on burn start and burn duration obtained in optimization 
process assuming impulsive burn are used to fix the bounds of the design variables in 
the optimization process with finite burn. These initial estimates help narrow down 
the bounds and hence reduce the computational time of the optimization process 
with finite burn. In the optimal finite burn strategy, the objective function includes 
minimizing the deviation in the intermediate orbit sizes, AOP during each burn, the 
TLI epoch and the final target orbit elements. In optimizing the above objective 
functions, suitable weights have been added to the terms of the objective function. 
Despite the smaller bounds, the computational time for the optimization process 
with finite burn is about half a day. In order to further reduce the computational time, 
multi-threading has also been implemented in the C++ code, developed to find the 
optimal burn strategy. 


4 Differential Evolution 


Differential evolution (DE) is a meta-heuristic procedure that optimizes a problem by 
maintaining a population of candidate solutions and creating new candidate solutions 
by combining existing ones using crossover probability and then keeping whichever 
candidate solution has the best score or fitness on the optimization problem at hand. 
DE can search very large spaces of candidate solutions, thus ensuring global solution. 
Three operations (i) mutation, (ii) crossover and (iii) selection are used to introduce 
variation in the population of candidates. The DE scheme is used in this study is 
rand/1/bin, which means that the vector to be mutated is selected at random, only one 
difference vector is used for mutation, and crossover is determined by independent 
binomial experiments. The performance of DE is influenced by three parameters (i) 
crossover probability (CR); (ii) population size (NP) and (iii) mutation factor (F). 
Further details of the DE algorithm can be found in [2]. 
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5 Algorithms 
5.1 Optimization Process with Impulsive Burn 


The unknowns in the impulsive burn strategy are the size of the intermediate phasing 
orbits around the Earth. The number of burns is fixed. 


1. Fix bounds on the sizes of the intermediate phasing orbits. If there are four 
burns, then there will four such bounds. 

2. Select values for sizes of the intermediate orbits randomly from their respective 
bounds. 

3. Propagate the orbit from injection epoch to an orbit pass which satisfies the 
visibility constraint from a network of tracking stations. Note that the prop- 
agation is not by using numerical integration. It is an approximate analytical 
process. After every orbital lap, Right Ascension of Ascending Node (RAAN) 
and Argument of Perigee (AOP) of the orbit are to be updated to account for 
J2 effect [1]. The other orbital elements are kept the same, and the state vector 
at perigee is computed. The check of visibility is done for this state vector. A 
constraint for an elevation angle of at least 70 at perigee has been imposed. This 
is done to later ensure visibility over the burn arc for finite burn. 

4. Introduce impulsive burns at the perigee to achieve the next immediate inter- 
mediate orbit size. The velocity requirement at the burn location is given 


by 
2 1 2 1 
Avg = ,/ Uh 
TB a rB a2 


where rg is the radius of the burn location and aj, a2 are the semi-major axes 
of the initial and target orbits. 

5. Repeat steps (3) and (4) till the number of burns (fixed aprior) is completed. 

Maximum number of orbital laps in the same orbit has been fixed to four. 

In case of non-visibility of spacecraft for four consecutive perigee of the same 
orbit, impose a heavy penalty in the objective function, which will eventually 
lead to rejection of the proposed solution. 

6. Find the apogee of the achieved orbit after all prefixed number of burns and 
compute IRachieved or Nargetl- 

7. Find the time elapsed to reach the perigee of the final intermediate orbit and 
compute ITLI Epochgchievea — TLI Epochiargedl. 

8. Check the visibility of the spacecraft at TLI epoch. In case of non-visibility, 
heavy penalty is to be imposed on the objective function. 

9. The objective function is set as 


F =w,X\|TLI Epoch 


+ Penalty for burn non-visibility at the perigees 


— TLI Epochyarget| + WX |Machieved = target | 


achieved 
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where w; are weights for the ith terms in the objective function. The weights in 
the objective function have been fixed after several trials and are very specific to 
a mission. For Chandrayaan I Mission later described in Sect. 6.1, the weights 
are w; = 500, w2 = 1, Penalty = 108. 

10. Build a population of a fixed size (NP) by repeating the steps (2) to (9). The 
population is a NP x (m + 1) matrix, where n is the number of unknowns, and 
the last column will be the value of the objective function. 

11. The DE optimization is carried out on the current generation, and the population 
is updated by repeating the steps (3) to (9) till the termination criteria are met. 


The output obtained from the impulsive burn strategy gives the size of the inter- 
mediate orbits, which in turn are used to find the burn start time assuming symmetric 
burn and duration of each burn for finite burn implementation. 

Given the initial mass of the spacecraft before each burn, ideal rocket equation 
gives the burn-out mass of the spacecraft after each burn: 


Minitiat 
Avg = gol in( ) 
sad M final 


where go = 9.81 m/s and J, is the specific impulse of the spacecraft. The burn time 
is given by 
__ Thrust 


n= 
8oLsp 


Minitial — Mina 
m 


Burn time = 


Thus, the estimated burn start time = Time of passage at perigee — estimated burn 
time/2, assuming symmetric burn. 


5.2. Optimization Process with Finite Burn 


The unknowns in the finite burn optimization process are the burn start time and 
burn duration of each burn. Previous Sect. 5.1 described how to find an initial esti- 
mate of burn start time and burn duration, which is used to fix bounds on these 
unknown parameters. The burn sequence and the intermediate orbit sizes are fixed 
to the values which are obtained in the optimization process with impulsive burn. 
These intermediate orbit sizes are used as targets to be achieved after each burn. 
However, achieving them exactly is not possible. During finite burn, the perigee 
altitude undergoes changes, which in turn changes the orbital period, and so, the 
target apogee altitudes after each burn also undergo small changes. The error in the 
achieved apogee altitudes can vary from 0.1% for the initial burns to 2% for the end 
burns. The objective function is formed to minimize these changes. 


1. Fix bounds on unknowns of all burns: burn start times and burn durations. 
2. Select values randomly from their respective bound for all unknowns. 
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Propagate the orbit from injection epoch to an orbit pass already identified with 
impulsive burn. Note that the propagation is using numerical integration of equa- 
tions of motion that includes perturbation due to J2 effect of non-spherical Earth. 
Unlike in the impulsive burn, the state vector gets updated in this numerical 
integration. 

Execute the finite burns at the selected burn start time for the selected burn 
duration for each burn. Check visibility of the spacecraft at all times during 
burn. In case of non-visibility of spacecraft during any burn, impose a heavy 
penalty in the objective function, which will eventually lead to rejection of the 
proposed solution. 

Find the apogee of the achieved intermediate orbit and compute INachievea,i — 
Ntarget,i| for the ith burn. Also, find the change in AOP of the orbit after the burn. 
Repeat the steps (3) to (5) for all burns. 

Find the true anomaly at the TLI epoch. 

Check the visibility of the spacecraft at TLI epoch. In case of non-visibility, 
impose heavy penalty on the objective function. 

The objective function for finite burn strategy is as follows: 


f = Wii (Error in AOP)) =P W12|Machieved, 1 ~ Narget,1| + Ww] (Error in AOP) 
+ W22|Machieved,2 — Ntarget,2| +...+ WLI (TAtLD 
+ Penalties for non visibilities during each burn and at TLI Epoch. 


where w;; represent the weights for error in AOP for ith burn and w;2 represent 
the weights for error in apogee altitude for ith burn. The weights in the objective 
function have been fixed after several trials and are very specific to a mission. 
For Chandrayaan I Mission later described in Sect. 6.1, the weights are w1, = 
105, wy2 = 2, w2, = 2000, w22 = 1073, w31 = 104, w32 = 10-4, wa, = 104, 
wri = 5000, Penalty = 1010. They can be easily set after studying the trend 
in the objective function on changing weights. 

Build a population of a fixed size (NP) by repeating the steps (2)-(9). The 
population is a NP x (7 + 1) matrix where n is the number of unknowns, and 
the last column will be the value of the objective function. 

The DE optimization is carried out on the current generation, and the population 
is updated by repeating the steps (3)—(9) till the termination criteria are met. 


6 Numerical Results 


A C++ code has been developed based on the procedures discussed in the sections 
above. To fix values for the DE parameters, a few trial runs have been carried out 
for an orbit-raising manoeuvre. The combination of the parameters which gave least 
error and minimal computational time is F = 0.4, CR = 0.5 and NP = 40. These 
set of parameters are specific to a sample orbital manoeuver, but they have given 
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satisfactory results for the other manoeuvers and missions discussed. To validate 
the formulation and the code, the following missions have been reproduced: (i) 
Chandrayaan I Mission and (ii) Mars Orbiter Mission (MOM). 


6.1 Chandrayaan I Mission 


India’s first Moon mission Chandrayaan-! carrying eleven scientific instruments for 
the purpose of expanding scientific knowledge about the Moon was launched on 
22 October 2008 from Satish Dhawan Space Centre, Sriharikota, India, by India’s 
PSLV-C11. The spacecraft was injected into transfer orbit of 254.4 x 22,932.7 km 
with inclination of 17.9° at 2008-10-22 01:10:19.081 UTC. The spacecraft was put 
into Moon’s polar circular orbit of about 100 x 100 km on 12 November 2008 by 
carrying out sequence of Earth and lunar bound manoeuvres. 

For this study, the spacecraft is assumed to be initially in EPO of 250 km x 
22,881 km. The following constraints were introduced in the optimization processes. 


1. The spacecraft should be visible during each burn from any of the ISRO ground 
stations, as mentioned in Table | (Fig. 1). 

2. No rotation of the apsis line of the intermediate orbits, during burn. 

3. Reach perigee of the last intermediate orbit at TLI epoch. 

4. Visibility at perigee of the last intermediate orbit for TLI burn. 


The parameters used in the study are given in Table 2. The optimization process 
assuming impulsive burn generated a burn sequence (4, 6, 7, 10, 11) which implies 
five burns at the 4th, 6th, 7th, 10th and 11th perigees of the orbit and the related orbit 
sizes. This sequence of burns and the intermediate orbit sizes are chosen ensuring 
visibility at the perigee of the orbit. The burn sequence and the intermediate orbit 
sizes are given in Table 3. The velocity impulse required and the burn duration 
corresponding to the intermediate orbit sizes and the mass of the spacecraft after 
each burn have been computed using equations mentioned in Sect. 5.1. The values 
have been tabulated in Table 3. 

It is to be noted this burn sequence was the same sequence proposed by ISRO. 
However, due to failure of burn at the 6th perigee, ISRO had to revise its burn strategy. 

The optimization process with finite burn is attempted using the results of impul- 
sive burn. The bounds on the start time and duration of each burn are fixed around 


eae ecu I Station Latitude Longitude 

Mission Trivandrum 8.536°N 76.78°E 
Bangalore 13.03°N T.51°E 
Biak 1.10°S 136.05°E 
Brunei 4.97°N 114.88°E 
Port Blair 11.64°N 92.71°E 
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Fig. 1 ISRO ground stations for tracking Chandrayaan I spacecraft 


Table 2. Chandrayaan I Mission specifications 


Specification Value Units 

Thrust 408 N 

Isp 318 

Initial mass 1380 kg 

PSLV injection epoch 22-10-2008 UTC 
01:10:19.081 

TLI epoch 03-11-2008 UTC 
23:15:00 

Target orbit 380,000 km 

Number of burns 5 


Table 3 Impulsive burn results for Chandrayaan I Mission 


Burn at perigee ha (km) Delta v (m/s) Mass (kg) Burn duration (s) 
4 37,393.45 318.37 1246.07 1023.69 

6 68,507.62 291.22 1126.35 831.70 

7 106,428.12 146.68 1070.49 388.14 

10 230,983.46 158.67 1013.17 398.22 

11 (TLD 380,000.00 57.11 993.29 138.06 

Total 972.05 993.29 
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the values, given in Table 3. The bounds on the burn duration are set within +40 s, 
whereas the bounds on the burn start time are different for each burn. It is set as 
+500 s for the first burn and increases to +1700 s for the last burn. The intermediate 
orbit sizes, along with few other parameters, have been used in the objective function 


of finite burn process, with suitable weights. 
To validate the formulation of the finite burn, the first burn of the mission has 
been reproduced and compared to the actual mission, as shown in Table 4. 
Small deviations in the current design from the Chandrayaan I Mission occur due 
to difference in the force model. The code accounts only for the J2 effect due to 
non-spherical Earth and not the forces, like solar radiation pressure, gravitational 
forces due to third bodies, other harmonics of Earth’s oblateness, etc. 
Table 5 gives the modified burn duration, the true anomaly of the burn start and 


Table 4 First burn design from this study and the Chandrayaan I Mission 


Orbit 3rd lap burn 4th perigee 4th lap burn Chandrayaan Units 
parameter start (with burn) end Mission 
hp 256.07 258.34 292.19 292.02 km 
ha 22,857.16 26,161.28 37,351.26 37,425.47 km 
a 179,314.61 19,587.81 25,199.73 25,236.75 km 
e 0.630097 0.661201 0.735307 0.735702 | — 
Inclination 17.89 17.89 17.89 17.89 2 
RAAN 353.09 35.08 353.06 353.07 e 
AOP 170.06 169.73 171.73 171.73 2 
TA 333.37 0.01 54.57 54.58 . 
Time passed 94,679.00 94,994.00 95,739.34 s 
since 
22/10/08 
01:10:19.081 
Burn 1060.34 1063.24 s 
duration 

Table 5 Finite burn strategy for Chandrayaan I Mission 
Burn at ha (km) | TA burn TA burn Delta v Mass (kg) | Burn 
perigee start (°) end (°) (m/s) duration (s) 
4 37,351 | 333.37 54.57 330.39 1241.27 1060.34 
6 68,507 | 323.05 33.87 287.16 1128.45 846.73 
vA 107,371 | 340.57 15.94 150.52 1075.28 397.43 
10 235,350 | 340.17 16.61 162.4 1020.72 404.73 
11 (TLD 380,000 | 351.98 4.27 54.6 1003 133.13 
Total 995.071 1003 
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Fig. 2 Visibility at different ISRO ground stations during different burns for tracking Chandrayaan 
I spacecraft 


the burn end of each burn. It also gives the modified orbit sizes and the corresponding 
velocity requirements. 

A single burn for TLI would require 1134.5 m/s against the 995.07 m/s with 
multiple burns. This could cut down the payload mass of the spacecraft by around 
40 kg. 

Figure 2 presents the plots for elevation of the spacecraft from different ISRO 
ground stations mentioned in Table 1. 


6.2. Mars Orbiter Mission (MOM) 


The Mars Orbiter Mission (MOM) is India’s first interplanetary mission. The Mars 
Orbiter Mission probe was lifted-off from Satish Dhawan Space Centre, using a PSLV 
rocket at 09:52.39.5 UTC on 5 November 2013. The MOM probe spent about a month 
in Earth orbit, where it made a series of six apogee-raising orbital manoeuvres before 
trans-Mars injection on 30 November 2013 (UTC) [5]. 

For this study, the spacecraft is assumed to be initially in EPO of 240 km x 
23,563 km. The following constraints were introduced in the optimization processes. 


Optimal Multiple Finite Burn Strategies for Trans-lunar ... 


Table 6 ISRO ground 


stations for MOM 
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Station Latitude Longitude 
Mauritius 20.35°S 57.55°E 
Bangalore 13.03°N TSE 
HBK 25.89°S 27.68°E 
Alcantara 2.40°N 44.41E 


1. The spacecraft should be visible during the beginning of each burn from any of 


the ISRO ground stations, as mentioned in Table 6 (Fig. 3). 


2. No rotation of the apsis line of the intermediate orbits, during burn. 


> 


Reach perigee of the last intermediate orbit at TMI Epoch. 


4. Visibility at perigee of the last intermediate orbit for TMI burn. 


The parameters used in the study are given in Table 7. The target is to reach 


Fig. 3 ISRO ground stations for tracking MOM Orbiter spacecraft 


Table 7 Mars Orbiter Mission specifications 


Specification Value Units 

Thrust 425 N 

Isp 315 s 

Initial mass 1340 kg 

PSLV injection epoch 05-11-2013 UTC 
09:52:39.5 

TMI epoch 30-11-2013 UTC 
19:30:00 

Number of burns 5 - 


358 G. Aggarwal and R. V. Ramanan 


Table 8 Impulsive burn results for MOM 


Perigee ha (km) Delta v (m/s) Mass (kg) Burn duration (s) 
0 23,563.37 1340.00 

5 28,224.61 123.41 1287.52 381.43 

8 40,201.14 214.31 1201.23 627.22 

10 72,605.17 272.71 1099.73 737.75 

13 119,733.40 155.77 1045.65 393.09 

15 190,533.00 97.27 1013.24 235.59 

Total 614.75 837.42 


an apogee altitude of 192,874 km from an initial apogee altitude of 23,563 km. 
The optimization process assuming impulsive burn generated a burn sequence (5, 
8, 10, 13, 15, 19) which implies six burns at the 5th, 8th, 10th, 13th, 15th and 19th 
perigees of the orbit. The burn at the 19th perigee would set the spacecraft in a 
hyperbolic orbit, along which the spacecraft escapes the Earth. This sequence of 
burns and the intermediate orbit sizes are chosen ensuring visibility at the perigee of 
the orbit. The burn sequence and the intermediate orbit sizes are given in Table 8. The 
velocity impulse required and the burn duration corresponding to the intermediate 
orbit sizes and the mass of the spacecraft after each burn have been computed using 
equations mentioned in Sect. 5.1. The values have been tabulated in Table 8. This 
results in velocity requirement of 1452.17 km/s, only for the impulsive orbit-raising 
manoeuvres to reach the target orbit before injection into trans-Mars trajectory. This 
burn strategy is the same as proposed by ISRO for MOM Mission. But ISRO could 
not implement this burn strategy due to incomplete 4th burn at the 13th perigee of 
the orbit. 

To find an optimal finite burn strategy, bounds on the start time and duration 
of each burn are fixed around the values, given in Table 8. The bounds on the burn 
duration are set within £150 s, whereas the bounds on the burn start time are different 
for each burn. It is set as +500 s for the first burn and increases to +20,000 s for 
the last burn. The intermediate orbit sizes, along with few other parameters, have 
been used in the objective function of finite burn process, with suitable weights. The 
optimal finite burn strategy obtained is given in Table 9. 

It is evident that the velocity impulse during finite burns for MOM is around 
904.62 km/s. 

The visibility plots, shown in Fig. 4, clearly show that the visibility constraint has 
not been rigidly been applied for MOM because of the rarity in the opportunity for 
the interplanetary mission. ISRO did eventually opt for partial blind burns to avoid 
delay or failure in mission. 
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Table 9 Finite burn strategy for Mars Orbiter Mission 


Burn at hq (km) TA burn TA burn Delta v Mass (kg) | Burn 
perigee start (°) end (°) (m/s) duration 
(s) 
5 28,237.38 | 341.02 13.81 124.11 1287.23 383.56 
8 40,501.48 | 330.32 25.085 220.31 1198.62 644.03 
10 73,750.47 | 325.65 29.91 278.42 1095.31 750.88 
13 120,222.31 | 340.29 14.97 151.59 1042.86 381.25 
15 236,694.48 | 344.98 13.64 130.19 999.824 | 312.82 
Total 904.62 999.824 


7 Multi-threading in Optimization 


Multi-threading is a specialized form of multitasking, which deals with the concur- 
rent execution of pieces of the same program. A multithreaded program contains two 
or more parts that can run concurrently. Multi-threading has been implemented in the 
scheme of differential evolution, only to evaluate the cost function of individuals in 
a specific generation simultaneously. Evaluation of the cost function is computation- 
ally dominating and thus demands parallelization. A sample lunar mission has been 
studied to highlight the importance of multi-threading to reduce computational time. 
The mission specifications are the same as shown in Table 2. The apogee altitude 
of the orbit is to be increased from 22,861 to 71,592 km in two burns. The mission, 
thus, has four unknowns—burn start time and burn end time of each burn. 

The computational time to find the optimal solution was reduced by eight times, 
after implementing multi-threading on a four-core processor (Table 10). The time 
taken to find the optimal five-burn strategy for Chandrayaan I Mission, with ten 
unknowns, took around 80 min with multi-threading. Without the visibility constraint 
for Mars Orbiter Mission, the computational time to find an optimal solution for ten 
unknowns is around 10—15 min. The computational time increases exponentially with 
increase in the number of unknowns. This brings out the need of multi-threading to 
reduce the computational time for an optimal solution. This computational time for 
the missions is aimed to be less than the orbital time periods of the spacecraft. This 
gives enough time to compute optimal finite burn strategy, in case of failure of a burn 
or any immediate necessary change in burn strategy. 
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Elevation v/s time for First Burn 


Elevation (degree) 


Elevation (degree) 


1.223 1.224 1.225 
Time passed since PSLV Inj. Epoch (secs) x 10° 


Elevation v/s time for Third Burn 


2.97 2.972 


Elevation (degree) 


Elevation (degree) 


Elevation (degree) 


2.974 2.976 
Time passed since PSLV Inj. Epoch (secs) x 10° 


Elevation v/s time for Fifth Burn 


0 
9.037 9.0375 9.038 9.0385 9.039 9.0395 
Time passed since PSLV Inj. Epoch (secs) x 10° 
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Elevation v/s time for Second Burn 


tt] 
2.102 2.103 2104 2105 2106 2.107 2.108 


Time passed since PSLV Inj. Epoch (secs) x 10° 


Elevation v/s time for Fourth Burn 


5.651 


5.652 5.653 


5.654 


Time passed since PSLV Inj. Epoch (secs) x 10° 


Port Blair 


Fig. 4 Visibility at different ISRO ground stations during different burns for tracking MOM 


spacecraft 


Table 10 Computational time for lunar strategy using single and multi-threading 


Unknown variables 


Single-thread 
Multi-thread 


1st burn 
start time 


(s) 
94,707.051 
94,707.051 


2nd burn 
start time 


168,384.810 
168,384.810 


Ist burn 
duration 


861.07 
861.07 


2nd 
burn 
duration 


966.61 
966.61 


Computational 
time (min) 


11.13 
1.34 


Generations 
in DE for 
computation 


178 
178 
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